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Abstract. For a variety of quantum dots, the widths of different single-particle levels 
may naturally differ by orders of magnitude. In particular, the width of one strongly 
coupled level may be larger than the spacing between other, very narrow, levels. We 
found that in this case many consecutive Coulomb blockade peaks are due to occupation 
of the same broad level. Between the peaks the electron jumps from this level to 
one of the narrow levels and the transmission through the dot at the next resonance 
essentially repeats that at the previous one. This offers a natural explanation of the 
salient features of the behavior of the transmission phase in an interferometer with a 
QD. The theory of this effect will be reviewed with special emphasis on the role of 
the interactions. New results on the dot-charging measurements and the fine structure 
of occupation switchings will be presented, accompanied by the unified description of 
the whole series of CB peaks caused by a single broad level. We then discuss the case 
where the system approaches the Kondo regime. 
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1. Introduction 

The challenging results of the experiment pQ which employed the double slit interference 
scheme in order to determine the phase [2] of the wave transmitted through a quantum 
dot (QD) in the Coulomb-blockade (CB) [3j regime is discussed in several articles of this 
volume. The unexpected result of this experiment was a repeated fast jump of the phase 
by — 7r between the resonances near the minimum of the transmitted current. [The phase 
at each CB peak increased smoothly by ir in accordance with the Breit-Wigner picture.] 
In this paper we will first review the mechanism intended to explain these experimental 
findings suggested in our paper Ref. [6] and then discuss new features and results. 
The essence of our suggestion in Ref. [6] was the observation that if a QD happens to 
support a single level extremely well coupled to the leads, the conductance for a series of 
charging resonances will be governed by the transport through (population of) this level. 
Other, narrow, levels are occupied, and the broad level is repeatedly depopulated via 
(almost) abrupt occupation jumps in the CB valley. This interlevel electron-exchanging 
mechanism, which has later been dubbed "population switching" (although the notion 
"level-occupation switching" may better reflect the physical situation), can lead to a 
satisfactory explanation of the experiment. The first suggestion of occupation switching 
as a possible explanation for the phase lapses and peak correlations (see below) was 
given in Ref. [7] and followed up in Ref. [E]. However, that mechanism was due to an 
assumed special structure of the dot, not to electron correlations as in the present work. 

Another important result of the experiment pQ, which also disagrees with a number 
of models, is the similar shape and height of consecutive CB peaks. This effect, which 
is often not emphasized enough, is again naturally explained within our picture. 

Experiments of the Weizmann group [H HI EJ dO, EH] measuring the transmission 
phase in QD-s have stimulated a wide theoretical interest, which is (partly) reflected by 

Rets. [TjEiiaiiaiiaiiiiisEaiH 

In particular, the phase behavior observed in Ref. p]: an increase by ir at the resonance 
accompanied by a sharp — ir jump in the CB valley, is in evident contradiction with what 
one would expect had the transmission of the current proceeded via consecutive levels 
in a 1-dimensional quantum well. In a two-dimensional QD the phase drops associated 
with the nodes of the transmission amplitude arise already within the single-particle 
picture [HI [151 US], related to the Fano effect. However, in order to have a sequence of 
such events one should consider a QD of a very special shape or level structure, which 
is not expected in the generic case [IB] . The model of Ref. [19] also does not allow 
to explain the series of drops, especially at very low temperatures. The mechanism of 
Refs. [74 [8] , as we already mentioned, makes nontrivial assumptions on the geometry of 
the QD and the way it changes under the change of plunger gate voltage. A generic 
mechanism suggested in Ref. [20] may indeed lead to the correlations in transmission at 
many consecutive valleys, but the predicted phase behavior differs from what has been 
seen experimentally. 

The mechanism for charging of the QD, which we will discuss in detail in the next 
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section includes two important ingredients. One is the necessity of the existence of levels 
anomalously coupled to the leads in sufficiently large irregular QD-s. The mechanism 
suggested for this purpose in our paper Ref. [6] will be reviewed in Sec. 2.2. Further 
developments of this and related mechanisms may be found in Refs. [13, 29, 32J. In 
all these cases the broadening of the single particle level does not require the electron 
interaction in the QD. 

The second ingredient is the interaction-induced occupation switching, which we 
describe in section 12.11 Further developments or applications of this effect, also 
discussing applications not related to the transmission phase problem, may be found 
in Refs. [2HI E21 EH ES EHJ EH EE] . Finally, recent papers [HI [IE] investigated, invoking 
the Fano effect and the dynamically generated level broadening [32J, both the universal 
regime, where phase lapses generically occur, and the mesoscopic regime, where the 
phase behavior is irregular, as in Ref. [TT] . The crossover between these two regimes as 
a function of a single parameter, was also discussed. 

Our main effort in the first parts of this paper will be to explain clearly the basic 
physics of the effect of level-occupation switching in the simplest cases, precisely defining 
the effect of the interaction. Effects of the electron spin [21] will be discussed in Sec. 3, 
while new results on the details of occupation switching will be treated in Sec. 4. The 
case where the QD approaches the Kondo regime [22] will be considered in Sec. 5 in 
relation with the experiment of Ref. [9] . We conclude with some remarks on where the 
theory stands vis a vis the explanation of the experiments. 



2. A model for the level-occupation switching 

A useful theoretical model for the description of charging effects in QD-s is the tunneling 
Hamiltonian in the constant interaction (Ucb) approximation (see e.g. [3]) 

H = ^ eicfci + U CB ^ (1) 

i i<j 

+ 5>(£0 L 4 + 4 + E^ C K + H - c -] + L " R ■ 

k k,j 

Here c(c + ) and a(a + ) are the annihilation(creation) operators for electrons in the dot 
and in the lead and Si,e(k) are the single-particle energies. We do not introduce the k 
dependence of the tunneling matrix elements tj. Summation over spin orientations is also 
easily included, but we first treat the spinless case (which is experimentally realizable 
by applying a strong enough magnetic field to the dot. To avoid the orbital effects of 
the field, it may be applied parallel to the plane of the dot). Under the assumption of 
capacitive coupling to the gate, the levels in the dot flow uniformly with the voltage 

Si = ei{V g = 0)-V g . (2) 
The energy of an electron in the wire with momentum k is given by 

*(*) = ~E F . (3) 
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Although the experiment [I] was clearly done in the CB regime, the widths of the 
resonances turned out to have been anomalously large, only few times smaller than the 
charging energy and larger than the dot's level spacing. Also, the widths and heights 
of all observed resonances are very similar. These surprising features of the results of 
Ref. P, which have not attracted so wide an attention as the phase jumps, motivate 
us to consider the model in which the QD supports a single very broad (well coupled 
to the leads) level, which we describe in the next subsection. Widths of all other levels 
will be assumed to be small and may be set to zero in the leading approximation. 

It is generally believed that the CB is observed only if the widths of resonances are 
small compared to the single-particle level spacing in the dot A [37]. This condition 
assumes that the couplings of all levels to the leads are of the same order of magnitude. 
However, as we will show below in section 2.2, for many situations the widths of the 
resonances may vary by orders of magnitude. In this case it does not make sense to 
compare the width of few broad resonances with the level spacing, determined by the 
majority of narrow, practically decoupled, levels. [Even the case of a single level being 
wider than the charging energy Uqb m &y be considered within our approach [33] .] 

2.1. The case of a single broad level 

Now we turn to the many-particle effects arising for the Hamiltonian (CQ) in the case of 
only one (iV-th) level in the dot coupled strongly to the wires, \t^J R \ 3> j ^ N. If 

the width of this level is larger than the single-particle level spacing A, a very nontrivial 
regime of charging of the QD may be described by means of second-order perturbation 
theory estimates (compare with Ref. [38]). This simple treatment produces all the 
important features of the problem (a similar approach was used recently in Ref. [39] for 
the calculation of CB peaks positions). 

The width of the well coupled level is given by 



The widths of the other levels are taken to be much smaller than the level spacing and 
may be neglected in the first approximation. The charging energy is still very large 
Ucb r. We shall show that transmission of a current at about (T / A)Iti(Ucb /T) 
consecutive CB peaks will proceed through one and the same level En- Let the levels 
with i < in the QD be occupied. Our aim is to find the total energy of the true 
ground state of the dot at different values of V g . Without loss of generality we may 
assume that the summation over % in Eq. (CQ) goes only over i > 0. (Thus we subtract 
from the total energy the trivial constant corresponding to selfinteraction of electrons 
with i < 0. Coulomb interaction between electrons at the levels with i < and i > is 
included into e^q.) We also subtract from the total energy the trivial energy of electron 
gas in the leads Yl l £ (k){ a k a k) ■ 

In this section we consider spinless electrons. We use a superscript (n) in the total 
energy of the many-electron ground state to denote the number of electrons, n, in the 




(4) 



i=L,R 



Level- occupation switching and phase anomalies... 



5 



LU 



CD 
O 
=) 



-0.2 - 



-0.5 - 



-0.4 - 



-0.3 - 



-0.1 r 



0.2 - 



0.1 - 



- 




-0.3 -0.2 -0.1 







0.1 0.2 0.3 0.4 0.5 




Figure 1. The evolution of the ground state energy of the QD with a single level 
strongly coupled to the leads, described in the text. Dashed lines show the energies for 
different states of the dot if the coupling to the leads would be switched off completely. 
The thin horizontal line is the energy of the empty dot, E = 0. Tilted dashed lines 
show the energy of the decoupled QD with the electron in state j, E = £j (lower 
line), and with the electron in state N, E = en (upper line). The ground state of 
the decoupled QD follows the thick dashed line. After the coupling of the level N 
to the leads is switched on the energies of the states having initially no electrons in 
the dot, E^l, and one electron in the state j, E^l, are renormalized, as described by 
Eqs. (|5l6l7l8p . The crossing of the two solid lines in the right lower part of the figure 
shows the occupation switching, Eq. ([9]). 



narrow levels (with i > 0) in the QD. Thus E^J and E\J. stand for zero and one electron 
in the narrow levels, respectively. We start with (0) and with the case of large positive 
£/v(^g) r. The only nontrivial contribution to the total energy is given by the second 
order correction where an electron hops from the lead to the level N and back, 



To calculate the integral here we use that (dn/de)de = (L/ir)dk, where L is the length 
of the wire and the energy and momentum are related by Eq. ([3]). In Eq. (jSJ) and 
everywhere below, £jv,£j are functions of ([2D- When with the increasing voltage V g 
the decreasing energy level crosses the region \e N \ ~ T the level N becomes occupied 
(occupation of the level j is energetically unfavorable, as we will see below, even though 
Ej < En)- 

The total energy of the system (with all the subtractions described above) may be 
calculated similarly to Eq. (jSJ) for en <S — r (the broad level being well below the Fermi 
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energy). The two contributions to E^J are now the energy of the electron in the dot En 
(which was taken where from the state e(kp) = in the lead (j3J)) and the second order 
correction describing the electron hopping from the level N to the leads. This gives 

E®=s K -L*(pc). (6) 



27T \\Sn 

Alternatively, Eqs. (jSJ) and (jfJJ) may be derived without explicit reference to many- 
body perturbation theory. To this end we first perform a unitary rotation of the 
lead states such that only a single effective lead defined by the annihilation operator 
dk = (t La k + t R(1 k)/ V\t L \ 2 + l^ R l 2 remains coupled to the dot. In the absence of the 
QD level the levels in this lead are quantized via k n — nir/L, n — 1, 2, • • • , oo. Adding 
a discrete QD level to this dense grid of the lead levels results in a slight shift of 
the new levels with respect to the old ones (see e.g. the book [40] ) . The levels with 
e(k) < En are shifted downwards, while the ones with e(k) > En are shifted upwards. 
The individual shift of any level is smaller than the level spacing in the lead (oc 1 /-£/), 
but the accumulated effect of the shifts of many (oc L) levels lead to a visible ~ T 
correction to the total energy. The sum of the downward shifts of all (full) levels with 
e(k) < in case of En > (en 3> T) gives Eq. ([51). The second term in Eq. (jBJ) (the 
case of En < 0) contains two contributions: a negative one due to the shifts of levels 
with e(k) < En, and a positive due to the shifts of the levels with en < s(k) < 0. A 
perturbative calculation of this correction becomes possible because of the cancellation 
of contributions from regions 2e n < e(k) < En and En < s(k) < 0. 

Finally, a precise treatment for a single state interacting with a continuum for 
spinless electrons yields [6]: 

-r 



P (o) 



+ — cot — , (7) 

7T 1 



which coincides with the Eqs. (|5|6l) at |ejy| ^> T and gives the correct interpolation 
between them. 

Let us now consider the branch where the level j is occupied (which we label by 
the superscript (1)). The energy of the electron in the QD is Ej. However, the energy 
of the empty level N is now increased by Ucb, due to the interaction. Adding one more 
electron via the hopping from the Fermi level in the lead, to that level, now costs 
En + Ucb- The ensuing reduction of the downward shift of the level E^j is of crucial 
importance. This is precisely where the effect of the interaction is felt. The analog of 
Eq. © for V g < Ucb now reads 

Without the interaction Ucb, the two second-order correction terms on the RHS's 
of Eqs. and (jSJ) are equal. This explains why these corrections to the total energy 
of the system do not affect the order of level populations in the noninteracting QD. For 
the case considered, Sj < En, level j will be the first to be filled when its energy crosses 
the Fermi level in the leads, £j(V g ) = 0, as is shown by a thick dashed line in Fig. 1. 
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With the interaction Ucbi after the electron is added to the dot the second order 
corrections in Eqs. ([6]) and (jSJ) start to depend on which of the orbitals N or j is occupied 
by this electron. In a large part of the CB valley the energy gain because of the increased 
second order correction overcomes the loss, e^ — Ej, caused by the raising the electron 
from the level j to the level N. The true ground state of the system here will be E^j. ([6]). 
The two functions E^(V g ) and E$(V g ) cross at 

UcB /qn. 

£n ~ exp{27r(£ iV -£ i )/r} + l ' [ ) 

and the ground state jumps onto the branch E^ t . Here the "active" electron inside 
the dot is transferred from the broad level N to the narrow one j. The energy of the 
current-transmitting virtual state N becomes positive again, e n + Uqb > 0. Thus, the 
phase of the transmission amplitude has returned to what it was before the process of 
filling of level N and the subsequent sharp jump into the state where level j is filled. It is 
the latter jump which provides the sharp drop by tt of the transmission phase, following 
its gradual increase by tc through the broad resonance. Fig. [1] illustrates the behaviors 
of the bare levels and the renormalized ones, including the crossing of the latter at the 
switching point. 

For the case considered so far, where the width of the narrow level vanishes, the 
switching transition is infinitely sharp. However when the coupling of level j is switched 
on, second-order processes involving jumps from level j to the lead and back from the 
lead to level N (or vice versa) will create a matrix element connecting these two states. 
This makes the crossing of the energies E^(V g ) and E^H(V g ), considered above, an 
"avoided" one, and provides a finite width for the switching transition. 

Even when the switching transition is broadened, the phase jump should still be 
abrupt (has zero width) in case of spinless electrons, as was proven in Ref. [2T] . Fast 
phase drops in the valley are associated with nodes of the amplitude, A ~ e— e , where e 
is close to the real axis (small Im^o)- For spinless electrons, the linear conductance (even 
through the interacting QD, since we are interested in the elastic processes) is determined 
by a single complex amplitude, G oc \A\ 2 , which is the amplitude of transmitted electron 
wave Ae lkx in e.g. the right lead for unit amplitude of the incoming wave, 1 x e lkx , in 
the left lead. [In case of spin the linear conductance would be given by a sum of the 
squared moduli of several complex amplitudes describing the scattering for different 
orientations of the electron spin and the spin of the QD, including spin-flip amplitudes, 
while the Aharonov-Bohm phase will be determined by the phase of the sum of non- 
spin-flip amplitudes, as discussed in Sec. 5 below]. The vanishing of the transmitted 
electron wave ip = Ae lkx at e = Eq means the existence of a solution of the Schrodinger 
equation with ip = in one lead. If the Hamiltonian is time-reversal symmetric, the 
complex conjugated wave function should also be a solution with the same boundary 
condition and e = ejj, hence Itueq = 0. The phase changes abruptly by it when the 
real energy Eq crosses the Fermi energy. [For the phase jumps associated with the exact 
vanishing of the conductance there is no way to distinguish between +tt and — tt jumps. 
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Figure 2. A sequence of charging resonances described (exactly for spinless electrons) 
by the theory of Ref. [BJ. Upper frame: the phase in units of 7r; lower frame: the 
transmission, or conductance in units of e 2 /h. Note that both the number of repeated 
resonances, where the occupation switching takes place in the CB valley, and the 
number of asymmetric peaks, where the narrow level is charged inside the broad 
resonance, are parametrically large, ~ V / A. [For this figure we chose Ucb/F = 4 
and U/A = 25. Increasing the ratio T/A would allow having a larger number of 
repeated peaks in the central part of the figure.] The conductance was calculated by 
cutting the usual Breit-Wigner resonances (two of them shown dashed) " with scissors" , 
and shifting them. 



This ambiguity is resolved, after the phase jump acquires a some width [22j [26], due to 
e.g. a finite temperature.] 

For the occupation switch taking place not too close to the charging resonances, 
detailed calculations of the conductance and phase behavior, may be done within the so 
called "elastic cotunnelling" approximation, as we shall describe in Sec. 4.2. Resulting 
from the destructive interference of quantum amplitudes, the exact vanishing of the 
transmission amplitude here has some similarity to the celebrated Fano zero. The 
broadening of the phase lapses in the case of broken time reversal symmetry was 
investigated in Ref. [23]. 
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The level-occupation switching guarantees the universal occurrence of the phase 
jump. Otherwise, it depends on the relative signs of the matrix elements connecting 
the states of the dot to the leads [26j E3 EH]. For example, the zeroes of the amplitude 
do not occur in the one dimensional case, but should be found in about a half of the 
spin zero valleys in a two-dimensional QD. While interesting effects can be found by 
adjusting these signs, we believe that a more universal mechanism is called for. 

When there are a number of narrow levels in the vicinity of the broad level, the above 
process can repeat many times. From Eq. (jHJ), we see that the number of narrow levels 
which switch occupations with the broad one is on the order of (~ (T/A)Iti(Ucb/T))- 
Thus the above number of consecutive resonances are due to the transition via one 
and the same level N. Obviously, the conductance and the phase shift (including the 
7r jumps) will be very similar for all these resonances, as long as they occur near the 
middle of the CB valley. When the occupation switch occurs close enough to a CB peak, 
the shape of that peak will also be deformed (see below). For other peaks, the phase 
will increase in a Breit-Wigner way by 7r around each resonance and sharply jump down 
by 7r at the switching point between the resonances. There is really no need for "further 
calculation" of this phase behavior. 

An example of such a series of consecutive resonances is shown in Fig. 2. We 
note the typical increasing asymmetry of the conductance peaks while going away from 
the center of the CB valley. A peak asymmetry with a qualitatively similar structure 
has in fact been observed by Lindemann et. al. [41] . It was interpreted as due to the 
conductance being determined by a special "bouncing ball" state having a width much 
larger than the others [29]. The mechanism invoked appears however to be different 
from the one described here. 

For the tails of series of peaks shown in Fig. 2, the occupation switch takes place 
within the charging resonance. With increasing relative peak number, \j — N\ ^> T/A, 
the usual narrow Breit-Wigner peaks in the conductance due to a population of level j 
should be formed here. The shapes of these peaks, as well as (possible) Fano-type 
conductance zeroes and 7r-phase drops will be strongly affected by the many-body 
orthogonality catastrophe and their investigation goes far beyond the scope of this paper. 

For electrons with spin, the Breit-Wigner-related formula (GO) does not work, but the 
perturbative formulas such as (jSJ) and (jEJ) may be still used away from the resonances. 
We will discuss the new effects caused by spin in section [3j Before that we explain in 
the next subsection why our model is likely to apply for a generic (smooth) QD with 
mixed classical dynamics. 

2.2. Justification of the model 

We now discuss the conditions for obtaining levels with strongly varying widths in 
the QD. 

An example of a system for which the widths differ drastically is the integrable 
QD El. However, it is hard to believe that the large (JV e ~ 100 4- 1000) QD may 
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Figure 3. The density of electrons in the dot at the resonances coupled to a single- 
channel wire(attached from the left). The line V — is shown dashed. Numbers 
correspond to the number of the level in the QD. About 95% of the norm of the wave 
function in the dot is shown. The "twin copies" (such as 100 and 112) of levels 94, 96, 
101, 103 are not shown. 

be even close to integrable. Nevertheless, at least in classical mechanics, a considerable 
gap is left between integrable and fully chaotic systems. Even in a nonintegrable dot 
two kinds of trajectories - quasi-periodic and chaotic - may coexist. In this case, in 
2-dimensions any trajectory (even a chaotic one!) does not cover all the phase space 
allowed by energy conservation. Consequently, the corresponding wave functions do not 
cover all the area of the QD. If such a regime is realized in QD-s, it easily explains why 
the widths of the resonances may vary by orders of magnitude. Moreover, many other 
features of such a QD may differ strongly from those of the chaotic QD |42j. An explicit 
numerical example, which supports the existence of such a regime is given below [6] . We 
require the QD not to be fully chaotic (neither do we require an integrable QD). It is 
not clear, how chaotic was the dot used in the experiment. However, the QD containing 
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~ 200 electrons was ~ 50 times smaller than the nominal elastic mean free path. Thus, 
disorder should not have been essential for the dynamics of the electrons. 

To illustrate the relevance of the model Eq. (CQ) with single strongly coupled level 
we performed numerical simulations for a model QD of a size I with a simple polynomial 
potential (a smooth QD coupled to two leads) 

V(x, „) (i _ + (y + 1 (l + 8 (f - i) ') . (10) 

Due to the strong mixing of the x and y coordinates the dot is expected to be 
nonintegrable, but, similarly to the experimental geometry [Tj, it is approximately 
symmetric. For simulations we considered the QD on the lattice and used / = 10 which 
was equivalent to 50 lattice spacings. The kinetic term is given by the standard nearest 
neighbor hopping. Below we present the results of calculations with the hopping matrix 
element r = 18 which corresponds to the dot with ~ 100 electrons or ~ 200 if the spin 
is included (similar numbers to those in the experiment). We have used the potential 
V of Eq. (|T0|) for < x < I. The lead formed by the potential V = 3y 2 was attached at 
x < and a hard wall was put at x — I. Within the energy interval 1.5 < e < 4.7 
exactly one mode may propagate along the lead. The analysis of solutions of the 
Schrodinger equation within this interval allowed us to find the positions and widths of 
quasi- stationary levels in the dot. As we expected, the widths fluctuate very strongly 
from level to level (by many orders of magnitude). In particular the widths of two levels 
#102 and #108 exceed sufficiently the level spacing T/A ~ 6 (the number of states 
doubled due to spin). The origin of the hierarchy of widths becomes clear from fig. 1, 
where we have plotted \4>\ 2 in the QD for (real) e at the top of corresponding resonances. 
The quantized version of different variants of the classical motion may be found in this 
figure. The narrowest level #103 corresponds to a short stable transverse periodic 
orbit. Other broader levels, such as #96, 106, may be considered as the projections 
of the invariant tori corresponding to quasi-periodic classical motion. This classical 
trajectory reaches the line V(x,y) = e only at few points. The candidates for chaotic 
classical motion (e.g. #110) also correspond to relatively broad resonances. Even in 
this case only a part of the QD is covered by the trajectory. For the most coupled levels 
#102 and #108 the area covered by the trajectory touches the lead by its corner. 

Moreover, two well-coupled trajectories contribute to the level #102. This is seen 
from the fig. 2 where we show also the \ip\ 2 at the left and right wings of this resonance. 
One contribution corresponds to the strongly coupled quasi-periodic trajectory (left), 
having the "turning point" V(x, y) = e just at the left contact. The other contribution 
comes from the true periodic trajectory (right). Two quantum states in the dot 
become mixed via interaction with the wire and form one broad (#102) and one almost 
decoupled (#104) resonance [43]. [This mechanism was also found and investigated 
in ref. [32] and the relationship to the Dicke effect [44] highlighted. It will be further 
discussed at the end of this subsection.] We have repeated the calculations several times 
for slightly different V and in a broad range of variation of the hopping. Typically, 
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Figure 4. Decomposition of the level 102 into parts corresponding to simple classical 
trajectories. Numbers are the energies for which the figures were done. 



we saw resonances of very different widths and the origin of the broadest peaks was 
explained by simple classical arguments. The explicit example also allows to clarify the 
range of validity of our method. In addition to the evident conditions A <C Tn <C Uqb, 
we should also have rv <C A\fW ei N e being the total number of electrons in the QD. 
For ~ Av^/Ve even the few broad levels start to overlap and the CB is expected to 
be eliminated. 

Taking into account the different sensitivities of longitudinal and transverse modes 
to the plunger [3 E] may allow to keep our broad level En even longer within the 
relevant strip of energy. This may provide an explanation for the even longer sequences 
of resonances accompanied by the — 7r jumps. 

In a more refined approach, adding new electrons into the QD should cause a change 
of the self-consistent potential V(x, y). The total energy of the dot and the wire will be 
lowered in the presence of strongly coupled levels. This may cause the potential of the 
QD to automatically adjust to allow such levels, which will support our explanation of 
the experiment of Ref. pQ. 

The Dicke-type mechanism for the dynamic generation of anomalously broad levels 
was used for many years in nuclear physics [43]. In QDs the attention to this mechanism 
was recently attracted by Ref. [32]. The latter work also generalized the treatment to 
several channels. When several levels are coupled to a continuum and their widths 
are larger than their spacings, the couplings among them mediated by the continuum 
become important (the same coupling was alluded to as causing the avoided crossing of 
levels 1 and in section 2.1). One now has to diagonalize the (nonhermitean, due to 
the finite lifetime of the levels) effective Hamiltonian matrix. In the simplest case, when 
all n 3> 1 initial levels are degenerate, one of the resulting levels (the " superradiant" 
one [44J) will acquire an 0(n) width, while the n — 1 remaining levels will become 
very narrow. This will still hold qualitatively when the level spacings do not vanish, 
but are smaller than the widths. Such a generation of one broad level was also found 
in a more general treatment, including the interactions and employing renormalization 
group (RG) approaches, in Refs. [T71 fT8] . 
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Figure 5. The energy dependence of the probability to find the electron in the QD for 
fixed incoming current. Numbers of some of the levels in the dot are shown explicitly. 
The lower curve is the result of subtraction of the narrow resonances, which allowed 
to resolve the strongly coupled levels #102 and #108. 

We notice however, that our numerical example [6], described above, showed that it 
may not be easy to employ the Dicke mechanism in its full strength in a QD described 
by a realistic 2-dimensional potential. We were able to observe a formation of the 
superradiant state from two close resonances (levels #102 and #104 in Figs. 3 and 5). 
However, to arrange the situation with several overlapping resonances, necessary for 
this mechanism, is difficult since the averaged level width in realistic QD appears to be 
always smaller than the level spacing. 

3. Effects of the electron spin 

The goal of this section is to consider the new features which appear if one goes beyond 
the simplest model of the occupation switching considered in Sec. 2. This section will 
be concerned with the effect of the electron spin. Since adding the spin leads to a 
significantly more complicated picture, we consider in subsections 3.1 and 3.2 the case 
of a QD with only two levels, one well coupled and one negligibly narrow. [These two 
subsections present the results found in Ref. [2T] .] The experimental tendency towards 
producing small few-electron QD-s also encourages one to investigate systems with small 
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number of levels. The case of many levels with spin will be discussed in Sec. 3.3. 



3.1. Two levels with spin, level- occupations and the ground state energy 

Charging effects with "spinful" electrons are described by the same tunnelling 
Hamiltonian Eq. (pQ), where now all summations include the summation over the spin 
index as well. We consider a QD with two levels, whose energies depends on a gate 
voltage V g and level splitting A as in Eq. (J2J) 

ei = -V g , e 2 = -V g + A. (11) 

Only one level, ex, is well coupled to the leads, having the width (jlj) 

r! = r>A , r<cw (12) 

The dot is charged by one electron at V g ~ 0, Ucb, 2Ucb and 3Ucb- 

Our first aim will be to find the ground state of the system flll|12p at different 
values of V g in the limit r 2 — > 0. Similarly to Sec. 2, let us denote by E^}, E^}, E^j the 
total energy of the lowest state of the QD interacting with the leads, with the narrow 
level populated by, respectively, 0, 1 and 2 electrons (more precisely, E^ t is defined as 
the eigenenergy of the Hamiltonian ([1]) minus the trivial energy of the electrons in the 
leads ^2 £ {k)( a t a k))- The functions E^ t (V g ) evolve smoothly with the gate voltage and 
the (averaged) occupation number of the broad level 1 also changes continuously. For 
example, the branch E^ t corresponds to an empty level 1 at V g < 0, singly occupied at 
< V g < Ucb and doubly occupied at Ucb < V g . For t%' R = the functions E^ cross 
at some values of V g , which lead to sharp changes of the ground state. 

With the use of perturbation theory in t^ R it is easy to find E^ t far from the 
charging peaks. Below the first resonance (at V g <C — r) the true ground state is 
evidently E^}. However, already here the virtual jumps of the electrons from the wire 
to the level 1 give rise to the correction 

.(0) _o V- f° l^l 2 dnl - " r - 



=L,R 



E° = 2 £ / ^~^de = — In ^ , (13) 
^— ' ' - e — e N de n \ ex 1 



Ep 



i?S = £2 --m 



r , / AEt 



7i \6i + U C b 



E ( £ t =2e 2 + U C B--hx 



r / 4E r 



7T V^l + 2f/cB. 

The extra factor 2 in E^j compared to Eq. dSJ) accounts for the spin. It is clear that 
Efol in this region lies significantly below E^j and E^ t . When the (increasing) voltage 
crosses the region \V\ ~ T, the dot is charged by the first electron. However, this 
electron may stay in the dot on the level 1 (described by E^j) or on the level 2 
Depending on what level is occupied the perturbation theory gives, in this range of V, 



y 
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(2) r / AE r 



71 \e 1 + 2U CBi 

The first logarithm in E^j accounts for the virtual jumps of an electron from the level 1 
in the dot to the wire. The other logarithms correspond to virtually adding the second 
electron to the dot (having e\ + Ucb or e\ + 2Ucb instead of S\ in the denominator). 
The two levels E^} and E^ t cross at a gate voltage (compare to Eq. Qj) 



V ^ V ' = exp{-2?A/r } + 1 • (15 » 



This result is valid for both signs of A. For T ^> |A|, Eq. (fT5l) reduces to £\ ~ —Ucb/2. 
At V g = V 1 the electron in the dot jumps from the broad level to the narrow one. 

In the second valley, Ucb < V < 2Ucb the dot is charged already by two electrons, 
which again may populate the two doubly degenerate levels in the dot in different ways. 
Here one finds, 

E^ t =2e 1 + U CB - V -lJ v ^—\ , (16) 

7T \\s 1 + Ucb\J 

r r / AE F \ , / 4E F 



E£l=e 1 + e 2 + U C B- — Un l , * , I + In 



(2) _ , rr 1 1. ^ F 



Er o > = 2e 2 + U CB --ln 



2tt [ \\ei + U CB \J \e 1 + 2U CB 

r / 4E f 



71 V £ l + 2U CB 

First of all, we see that just after the charging peak, at V g > Ucb, the true ground state 
is Efg} . This is in contrast with the situation at V g < Ucb, where the ground state was 
EfH ( JT4l) . Thus we may conclude, that within the resonance not only is one electron 
gradually transmitted from the wire to the level 1 in the dot, but also a second electron 
is taken from the narrow level 2 to the broad one 1 (see fig. 6). In the limit r2 — > 
this second "transfer" is abrupt and takes place at some V g = W 1 ~ Ucb (f° r zero 
temperature, T = 0). This prediction of the possibility to have sharp features within 
the charging resonance is the main new effect caused by the spin. All the three energies 
(fIB|l cross at the same value of gate voltage (c.f. (TT5T) ) 

""^^^I^fiTT (17) 

and the true ground state becomes E\ [. Now already two electrons jump together from 
the broad level to the narrow one. At V g > V 11 ( TT5"j) the occupation of the quantum 
dot proceeds in a fashion symmetric to the above. At the third peak the third electron 



n(2) ^ 

^(2) 77(1) 



is added to the dot. Were the branch E^J the stable one, this would have been the 



uncoupled electron at the level 1. However, ElJ and ElJ cross at the top of the third 
peak (at V g = W 11 = W 1 + Ucb) and the ground state for the first half of last valley 
has a single unpaired electron at the narrow level. Finally, at 



= ^ + < 18) 
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Figure 6. Schematic drawing of the averaged occupation number (nj + n^) for the 
broad level 1 (solid) and narrow level 2 (dashed). The gate voltage V is measured in 
units of the charging energy Ucb- The charging resonances are at V g /UcB ~ 0, 1, 2, 3. 

E^ t and E^l cross and the broad level becomes singly occupied again. The fourth peak 
(V g ~ SUcb) completes the charging of two levels in the quantum dot by four electrons. 

To summarize the discussion of this subsection, we show schematically in Fig. [6] the 
averaged occupation numbers of our two levels (n^ + r^) and (n 2 j + n 2 [). The four 
charging resonances correspond to V 9 /Ucb ~ 0, 1, 2, 3. One may see on the figure many 
abrupt changes in the occupation of both broad and narrow levels. The occupation 
number of a given level is not measured directly in the experiment. However, the 
appearance of an occupation switching may be seen in the signal of a suitable QPC 
detector [15] as we describe in Sec. 14.11 In the following subsections we discuss the 
transport properties of the QD with a single broad level for "spinful" electrons. 

3.2. Two levels with spin. Conductance 

The zero bias conductance G of our two-level quantum dot is shown schematically 
in Fig. [7J In the limit of an "invisible" level 2 (T^ — > 0), we may introduce three 
conductances G^ 0,1 ' 2 ) corresponding to empty, singly-occupied and doubly-occupied 
narrow level (characterizing the three "ground state" energies E^l of the previous 
section). The role of the electrons in the level 2 reduces in this case to simply raising 
the current-transmitting level 1 via the Coulomb repulsion. Thus 

G^(V g ) = G^(V g + Ucb) = G^(V g + 2U CB ) • (19) 

We have shown schematically the function G^°> on the same Fig. [7J (vertically offset). 
In the limit T 2 Ti the curve for a two-level dot is obtained simply by cutting and 
horizontally shifting the parts of the curve for the single-level dot, in agreement with 
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Figure 7. The gate voltage dependence (schematic) of the conductance G. The upper 
curve shows G^ in the case of only one (broad) level in the quantum dot (vertically 
offset for clarity). The lower curve depicts the conductance given by the different 
in various regimes. Sharp features at the two peaks and three valleys are seen. 



Eq. (1191) . as shown in the figure. Thus the relations Eq. (1191) allow one to describe the 
singular behavior of the conductance even without an explicit calculation of G^°\ 

The sharp features seen in Fig. [7] will be smeared at finite temperature (see 
examples in Ref. [21]). To resolve them one needs to measure the conductance at 
T <C T. If the temperature is still large compared to the Kondo temperature, 
Tr = (UcbT /2y/ 2 e n£( - £+UcB ^ 2rUcB , the conductance away from the resonances may 
be found by calculating the transmission tunneling amplitude in the second order of 
perturbation theory ("elastic cotunnelling" ) . This gives, 

e 2 T L T R 



G (o) 



h 



e-2 
-I 



at V g < 0, 



(20) 



G (0) = -rT L T R 
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G (0) 



rvr 



£1 R 



U CB , 

0<Vg< U CB , 

at U C b < V g . 



(ei + u CB y 



h (s. + Ucb) 2 ' 

The first line in Eq. ( 1201) is for an empty dot. The second line describes the current 
through the QD having one electron on the level E\. The current in this case is 
a sum of three contributions [21] accounting respectively for the spin-flip processes, 
transmission of the electron with the spin parallel to the spin of the electron in the 
QD and transmission of the electron with antiparallel spin. Note that the spin-flip 
contribution vanishes at Ucb = 0, as it should. The last line in Eq. (|20|) is for the case 
of two electrons on the level E\. 
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Two kinds of sharp features are seen in Fig. [7J First, both a cusp and a jump take 
place at the peaks: W 1 , where the curve is replaced by and W n , where G^ 2 -* 
is replaced by G^\ The jump vanishes for A < T, but the pronounced cusp survives 
even for A = 0. Currently, no theory exists which would describe the fate of these 
singularities at T —>■ and finite (small) coupling to the leads of the level 62- 

Besides the above, there are three jumps of G(V g ) in the valleys. The values of the 
conductance at the borders of the first jump, at V g = V 1 ± are (assuming A <C T) 



The contribution from the spin flip processes for V g < V 1 is twice that from the elastic 
processes. Therefore, the conductance drops near V g = V I by& factor of 3. At V g = V n 
the two electrons jump from the broad to narrow level. The discontinuity at V 11 , 
which follows from the small difference in the probability of the electron-like and hole- 
like processes, vanishes for A <C T (restoration of particle- hole symmetry at A = 0). 
Examples which find the fine structure of such jumps at T 2 7^ will be given below. 

3.3. Many levels with spin 

A straightforward generalization of the theory presented in previous two subsections 
allows one to find the conductance (and level occupancies) for a multilevel QD with a 
single broad level, T ^> A, for "spinful" electrons. The resulting conductance behavior 
is shown in Fig. [8] (compare with Fig. 2). Similarly to the case of 2-level dot, considered 
in Sec. 13.21 the conductance curve may be found by cutting and gluing the parts of 
the curve G^{V g ) describing a QD with single (spinful) level. In contract to Fig. [JJwe 
now consider the case of zero temperature (T <C Tk). Thus the "generating function" 
for conductance of multilevel QD is now the conductance at T = for the Anderson 
impurity model [46], shown by thin dashed line in the figure. [An exact Bethe ansatz 
solution for this model [47] gives G^(V g ) explicitly] 

The series of peaks shown in Fig. [8] contains three large, ~ T/A, sets with 
qualitatively different charging scenarios. First comes a sequence of peaks (V g /UcB ~ 
0,1,2,3,4 in the figure) having one occupation switch per charging event. At the 
left wing of each resonance here the broad level in the QD is empty and is gradually 
populated with the increasing V g . Somewhere inside the charging resonance it becomes 
energetically favorable to populate one of the narrow levels in the dot by single electron. 
The broad level is depopulated at this moment and its energy is raised, En — > £n + Ucb- 
A reliable description of the fine structure of such an occupation switching for a finite 
coupling of the narrow levels to the leads, remains a challenging theoretical problem 
(similarly to the case of spinless electrons of Sec. 2.1). 

With further increase of the gate voltage the occupation switching, for a series 
of ~ rin(f/cs/r) peaks, splits into two parts. Now inside each peak one electron 
jumps from the narrow to the broad level in the dot. Such events for 2-level case were 
considered in previous subsections (V g ~ W 1 and V g ~ W n ). Thus within a single 
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and G « 16 
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Figure 8. A series of conductance resonances in a QD with a single broad level with 
spin. Here in the first five peaks electrons first try to populate the broad level with 
either spin up or spin down and then (before the charge at the broad level reaches 
1 x e) arc thrown away after the charging of the narrow level. Similarly the last three 
peaks (we don't show more because the picture is symmetric) the lead-electron jumps 
onto the narrow level in the QD, which raises the energy of occupied broad level and 
increases the transport through it. In the three peaks in the middle, one electron jumps 
from the narrow to the broad level in the QD inside the charging resonance. The jumps 
in the tails of the peaks here correspond to depopulation of both spin states of the 
broad level. The thin dashed line shows the (schematic) conductance for the Anderson 
impurity model at zero temperature, including the increase of the conductance upon 
adding the first electron, the Kondo plateau and the conductance drop after adding 
the second electron. The conductance for the multilevel QD is obtained by gluing and 
repeating the pieces of this curve. 



charging resonance the well-coupled level in the dot becomes occupied by two electrons. 
Between the resonances both electrons from the broad level jump to narrow ones. This 
is where one will see the n jumps of the transmission phase in the CB valley. 

Finally, for a series of ~ T/A peaks the charging events start (left wing of the 
resonance) in the situation when the broad level stays inside the QD (e^ ~ —Ucb) 
being populated by two electrons. Then, at a certain voltage, an electron from a lead 
jumps onto the QD narrow level. Repulsion from this electron raises the broad level, 
£n — * £n + Ucb ~ I\ and partially depopulates it. The occupation of the level N is 
completed again at the right wing of the resonance. 

The important consequence of the scenario presented in this section is the absence 
of the Kondo effect for a large series of charging resonances, which might be easy to 
verify experimentally. 



Level- occupation switching and phase anomalies. 



20 



4. Sharp features and charge detection 

In subsections 4.1 and 4.2 we considered the effects which appear if one takes into 
account the finite width of the narrow level. Since the particular features considered 
here do not depend on the spin the discussion is narrowed to spinless electrons only. 



4-1. Sharp features: QPC detector signal 

In this subsection we investigate the fate of the sharp features in the CB valley considered 
before, if the narrow levels in the QD acquire a small but finite widths. Let, as in 
Sec. 3.1, 1 be a broad level and 2 be a narrow one in the QD. First, in a way analogous 
to Eqs. ( l4lfT2~l) we introduce the width of the second level T 2 and the interlevel coupling 
Ti2 = 27r J^ i=L R t\t 2 *dn l jde. Consider a CB valley where there is one electron in the 
dot either on the level 1, or on the level 2. The dynamics of QD is described by the 
effective single-particle Hamiltonian, accounting for the coupling to the leads, whose 
matrix elements are [e 2 — £i = A) [21] 

„ T x (e + Ucb\ „ T 2 fe + UcB 

Hn = si - — In — j-j — , H 22 = e 2 - — In 



2tt — \ lei ) ' 2tt 



Hn= - T ^ln( £ -±^)=H; 1 . (22) 



2tt V \ £ \ 

Here Hn and H 22 are nothing more than the renormalized single particle energies E\ 
and e 2 , similar to Eqs. (I6HB1) minus a proper constant. In the argument of the logarithm 
e may either be E\ or e 2 , since e + Ucb, \ e \ 3> A. For simplicity we consider the time 
reversal symmetric case, H^ = H 2 \- Eq. (1221) is valid for an arbitrary ratio of couplings 
of two levels to the leads. In the general case however, the width of the level switch may 
become of order Uqb ( see Eq. (1281) and remark below it). Below we consider the sharp 
switch case, Ti > r 2 , T u ~ \/T\T 2 ~. 

For electrons with spin there are two species of the effective Hamiltonian Eq. (1221) 
for each spin orientation. [If in addition the QD has one electron for each spin, there 
are no spin-flip transitions in the linear conductance.] In this and in the next section 
we consider spinless electrons. 

Examples of resolving the sharp features in the conductance based on the avoided 
level crossing described by Eq. (1221) may be found in Ref. [21]. Here we give only two 
examples, which seemingly have not been covered in the existing literature. 

Using a side quantum point contact (QPC) as a charge detector [45] becomes 
now an efficient experimental tool for finding the number of electrons in small QDs, 
complimentary to the transport measurements. The same device may obviously be used 
to detect the occupation switching [36] . 

The ground state of the system described by the Hamiltonian (l22~l) is a superposition 
of two states 



U) = a |l> + m , *,„ = ±( /I T (23) 
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Figure 9. The signal at the detector (QPC) for two values of the parameters used 
in Eq. (j2"8|) . The dashed line shows the case of vanishing cross-scattering amplitude 
(B = 0). The solid line shows the case then the scattering amplitude for the QPC 
electron with flipping the state of the electron in the QD from 1 to 2, equals to the 
difference of direct scattering amplitudes for the dot electron in states 1 and 2{B = A). 



For weak coupling between the QD and QPC the transmission amplitude through the 
latter reads (St <C t ) 

t = r + St , St = (a 2 - ^ (W) ~ (^l 2 ) + 2a0(l\R\2), (24) 

where R is some single-particle operator and both the background transmission 
amplitude t and the matrix elements of R weakly depend on the gate voltage V g . 
Thus we find the change in the QPC conductance 

SG D = A - H ^~ H ™ +B - 2Hl2 , (25) 

v /(iJ 11 -^22) 2 + 4|i/ 12 | 2 ^/(H u -H 22 )i + 4\H 12 f 

where A = ^t ((1\R\1) - (2\R\2)) and B = ^2r (l|i?|2). To see explicitly the gate 
voltage dependence of SGd we consider the case of a single broad level, Ti 3> r i2 ^> T2- 
Now we may omit T 2 in Eq. ( |22l) and use an expansion around the center of the switch 
(ei = -V g ) 

£ ° - exp(-2 7 rA/r 1 ) + l' [ } 

which allows us to write 

Hn - H 22 w ^— sinh ( 2n^-) (e 1 - e Q ) , H 12 « ^A. (27) 
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Now, instead of Eq. (1251) one has 



+ 



kB 



&r 12 u C B 



(28) 



y/fa ~ E ) 2 + K 2 



y/(Ei ~ S ) 2 + K 2 ' 



r? sinh ( 2vrA 
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The width of the occupation switch, k, is always a fraction of the charging energy, 
becoming k = Vl2 ^ B for A << r\. The function 5Gd(£i) is shown in Fig. Note the 
asymmetry of the step in the QPC signal and the maximum on one side of the step for 
B 0. 

4-2. Sharp features: Phase drops 

While we declared in Sec. 2.1 that the phase lapses in the CB valley for spinless electrons 
should be abrupt and equal ir, the actual phase drops in Fig. 12.11 are always somewhat 
smaller than n. In this section we show how this exact 7r is restored after the coupling 
of the narrow levels is switched on (still we assume Ti r 12 T 2 ). We also continue 
working with spinless electrons here. 

To find the conductance and phase behavior we first notice that after having 
successfully employed the coupling of the narrow level to smoothen the occupation 
switch, there is no need to consider the contribution to the conductance from this small 
coupling. What we have now is a QD with two levels |1) and |2) ( 1231) . each coupled to 
the leads via tunneling matrix elements at L,R and (3t L,R . One level, |1), is occupied by 
the electron (energy e < 0), the other level is raised due to electron repulsion (e + Ucb)- 

In the valley between two CB peaks, the conduction is dominated by the so-called 
elastic cotunneling (i.e. second order) processes through the resonances at larger and 
smaller energies. For example, tunneling first from the left lead to the dot (whose 
relevant state must then be empty beforehand) then from the dot to the right lead. 
This "electron process" is dominant near the right-hand CB peak. The other possibility, 
dominant near the left-hand CB peak, where the relevant dot level is full, involves first 
a transition of the electron from the dot to the right lead, then from the left lead to 
the dot. Therefore, the latter process is termed a "hole process". For fermions, the 
hole process is well-known to acquire an extra minus sign. In our case, where due 
to the level-occupation switching, the two CB peaks are due to the same level, the 
transmission amplitude must have opposite signs in the two sides of the Coulomb valley. 
Thus the transmission amplitude in the time-reversal symmetric case must have a zero 
somewhere in-between. A simple calculation with the use of Eq. ( 1231) now gives (we 
assume t L = ±t R ). 



To extract the transmission phase from this result we further notice that since the 
couplings of two levels to the leads are proportional we may perform a standard unitary 
rotation of the leads [H], oc t L a k + t R a R , a' k oc —t R a k + t L a R , to define one coupled 
(a k ) and one decoupled (a' k ) lead. All transport properties are found from the S'-matrix 




(29) 
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Figure 10. The extended part of Fig. 2 showing the fine structure of the conductance 
and the phase behavior at the occupation switching, as described by Eqs. (|2U)) and 
(150)) . The phase jump now becomes exactly A0 = tt, while the conductance is smooth 
and vanishes at a certain point. The value of n (|28p is taken to be k = [/cs/50. 
The dashed lines show the bare curves obtained "with scissors", as in Fig. 2 (with the 
abrupt jump in the conductance and the phase jump A</> < n). 



for the single coupled lead, 5* = e 2l< ^. In particular the transmission amplitude and 
conductance are (again \t L \ = \t R \) 

e 2 

t = sin0e^ , G = — sin 2 0. (30) 
h 

This allows us to relate the conductance (l29l) and the phase of the transmitted wave <fi. 
Fig. 14.21 shows the conductance and the phase behavior found from Eqs. f l29l) and ( 1301) . 

An elastic cotunnelling treatment of transmission zeroes in noninteracting QDs 
was presented in Ref. [26J. Ref. [24J investigated the evolution of phase lapses in an 
interacting QD including the case of broken time-reversal symmetry. 

5. The sensitivity of the transmission phase to Kondo correlations 

Among the observations [19J [50j EH [52j [53] of the Kondo effect [M] in quantum 
dots(QD) [3B1 E2, EE], two experiments 0, ED] were devoted to the measurement of 
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the phase shift of the transmitted electron. These experiments were aimed at a direct 
observation of the fundamental prediction [57J EE] of the Kondo model: the <p = n/2 
phase shift experienced by the scattered electron after the spin of the impurity is screened 
into a singlet. 

In agreement with the general predictions and the numerical renormalization group 
calculations (30] for an Anderson impurity [16] in Aharonov-Bohm(AB) interferometer, 
the development of a plateau of <p i n the Kondo valley was seen in Ref. [9] (although 
the reported saturation value of the phase was not vr/2). A very important, unexpected 
feature of the experimental results was the strong sensitivity of the phase to Kondo 
correlations: the phase saturated at <fi ~ const, while the conductance in the valley was 
still well below the unitary limit, indicating the absence of Kondo screening. 

Our aim in this section will be to explain this strong Kondo effect in the phase 
[6]. We will see that although the nontrivial phase behavior is indeed governed by the 
Kondo physics, the main changes of phase take place at temperatures parametrically 
larger than T#. In a sense, we show that the phase changes in the regime where, 
although the spin of the impurity is not screened, the running Kondo coupling constant 
exceeds parametrically its bare value. 

The transmission phase is measured, as before, by embedding the quantum dot into 
one arm of the two- wave AB interferometer [1] . Let A d s be the transmission amplitude 
for an electron with spin projection s through the QD having a spin projection 5*. 
Respectively, let A r s be the transmission amplitude through the second reference arm. 
(Obviously, the spin- flip processes occurring in the QD and not in the reference arm 
do not contribute to the interference.) Now, the part of the current oscillating with the 
change of the magnetic flux threading the interferometer takes the form 

G AB *ReJ2 A r s *A d sS = ReA r * £ A d sS . (31) 

sS sS 

We use the fact that (for the weak magnetic field used in the experiment) the 
transmission through the reference arm does not depend on spin. By measuring the 
relative phase of the AB oscillations at different values of the gate voltage one extracts 
the information about the transmission phase. We see that in fact in the interference 
experiment not a single amplitude is measured, but the sum of all possible amplitudes 
corresponding to different orientations of the spins of both the transmitted electron and 
the QD [19]. 

5.1. The single-level Anderson model. 

Far from the charging resonances the interaction of the lead electrons with the dot is 
described by the Kondo Hamiltonian 

Hk = Y1 £ * c ?A + S d + VS*n«] . (32) 

iks ij 

Here i and j denote left(L) and right (R) leads, the operator c\ s annihilates the electron 
in lead i with momentum k and spin s. The Pauli operators and the density of the 
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conduction-electrons on the impurity are given by 

ff *e = C ks*ss<4>s> > = Yl C ks4's ■ ( 33 ) 

kk'ss' kk's 

Explicit formulas for Jo and Vq will be found below from the tunnelling Hamiltonian 
describing the QD. We will consider a time-reversal symmetric system. Then the 
matrices J % \ V %: > are real and symmetric. 

The interaction of the conduction electron with the spin of the dot may be 
diagonalized by an orthogonal transformation [59] of cl,cr into new operators c u and 
c v , described by the angle 9, tan26> = 2J LR /(J LL - Jq R ). This gives 

Has = {Jftu + J%K)S d . (34) 
As long as Jq , Jq -C 1 the two couplings are renormalized independently 

l fju,v = yjgv + Ap j n _ ^ 

Here v = = vr is the density of states in the leads. In the case of antiferromagnetic 
coupling (J > ) this formula may be written in the form 4z/J = l/ln(T/7V) with 
Tk being the Kondo temperature. Crucial for the understanding of the phase behavior 
is that in the leading order only the spin-dependent part of the Hamiltonian fl32l) is 
renormalized, while the scalar coupling remains unchanged, V = Vq = const. 

In the simplest case of only one level in the dot (Anderson impurity model), only 
one mode c = {tici + tRCR)/t is coupled to the dot via the tunnelling matrix element 
t = a/I^lI 2 + |^i?| 2 ) while the second mode remains completely decoupled. The bare 
values of the coupling constants in the Kondo Hamiltonian (132]) are now given by the 
second order of perturbation theory 

^ H + „/. r = Vo-#- (36) 



-e d 2 ' -(U CB + e d ) 2 

Here e d < is the energy of the impurity level and Ucb is the charging energy and the 
Fermi energy is Ep = 0. 

The non-spin-flip transmission amplitudes for parallel and antiparallel spins of the 
dot and the electron are 



-e d \ 
.4", -— ( ! 



^ = ^1-^), (37) 



U -(UcB+e d )\ -2(U C B + e d ), 

where T = 2nt 2 v. The imaginary part of the amplitudes here is formally of the fourth 
order in the tunnelling amplitudes tL,n- However, the calculation of the phase of 
transmission amplitude requires only the S-matrix for non-spin-flip scattering to second 
order in t. That is why the Kondo correlations do not contribute to the phase in 
the leading order and Eq. ( |371) coincides with the expansion of the usual Breit-Wigner 
formula. The Kondo effect appears in the real part of the amplitudes at the order ~ t A , 
which we will take into account via the renormalization fl35|) . In order to find the AB 
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Figure 11. The phase as a function of the depth of the impurity level — £d for 
r = U/m. Solid lines show the calculated <p(e d ) for T = T, T/4, T/4 6 . Since the 
theory is valid only for T ^> Tk, the curves are shown only for T > 2Tx(ed)- The 
dotted line depicts schematically the expected 4>( £ d) for T = T/1024. The thin dashed 
line shows the phase for the sum of two Breit-Wigner resonances and the thick dashed 
line is the Bethe-ansatz solution for T = 0. 



current one should rewrite the sum of the two amplitudes and in terms of scalar 
and magnetic couplings ( |36l) and then replace Jo by J ( 1351) . which gives 

G AB oc ReA r *[V Q - mv(V£ + J 2 /4)] . (38) 

This formula is the central result of this section. It should be compared with a usual 
conductance of a Kondo quantum dot, obtained by adding |A^| 2 , \A^\ 2 ( |37|) and the 
corresponding spin-flip contribution, 

G = 2^r 2 M^ (V 2 + \.J^j . (39) 

Remember that Vo is not renormalized, while J ~ 1/ \n(Y /Tk)- Eqs. ( |39j) and ( |38l) are 
justified only for vJ <C 1. This does not allow for the quantitative description of the 
conductance in the unitary limit, where G ~ 2e 2 /h. On the other hand, the phase shift 
close to 7r/2 develops at 

uJ < vJ ~ \/uVo < 1- (40) 

This gives the temperature scale, explaining the high sensitivity of the transmission 
phase to the Kondo effect observed in the experiment |9j: In (T/T K ) ~ ^J\n(T /Tk)- 

Fig. 11 shows the phase 4>{ed) found via Eq. ( 1381) for Y = Ucb/30 and different 
temperatures. Due to a node of V at e<i = —U/2 the phase equals exactly 7r/2 in the 
middle of the valley. The width of the phase drop at higher temperatures (T ~ T) 
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is ~ T. The exact Bethe-ansatz solution [UJ for T -C 7# predicts a broad plateau 
between resonances with <fi « 7r/2 [30]. However, due to a strong dependence of T K 
on the position of the impurity level Ed the phase changes from = O,7rto0 = 7r/2 
in a rather nonuniform way. For intermediate temperatures the phase first develops a 
shoulder with (f> « 7r/2 (T <C T# (£,$)), while in the center (T ^> Tx(^d)) pronounced 
maximum and minimum are formed to the left and right of the point Ed = —Ucb/^- 
This structure may be seen on the experimental figures of ref. [9], but appears not to 
have been discussed in the numerical calculation of ref. [30J. The conductance ( 139]) in 
the middle of the valley for the same parameters of the dot as in Fig. 1 varies from 
G = 0.013e 2 //i for J = J to G = 0A2e 2 /h for the lowest temperature T = T/4 6 . We 
see that the phase behavior is well developed when the conductance is sufficiently below 
the unitary limit. 

6. Conclusions 

To conclude, we have considered the model for which upon increasing the plunger gate 
voltage V g , it is energetically favorable to first populate in the dot the level strongly 
coupled to the leads. At a somewhat larger V g a sharp jump occurs to a state where 
the "next in line" narrow level 1 becomes populated. This jump occurs due to the 
appearance of Uqb m some of the denominators of the second-order energy shifts. It 
is thus a genuine correlation effect and can not be faithfully addressed via a Hartree- 
Fock-type approximation. The jump accounts for the sharp decrease by ~ 7r of the 
transmission phase. The similar strengths of resonances seen in the experiment [1] and 
their large width are also clear within our mechanism. It appears that this is the only 
available model which can explain these two important features of the experiment. Even 
independently of that, we believe that the level-occupation switching is interesting in its 
own right. It can be further studied experimentally using a combination of charge sensing 
and transport techniques. The current transmission through such a QD resembles the 
behavior of rare earth elements, whose chemical properties are determined not by the 
electrons with highest energy, but by the "strongly coupled" valence electrons. The 
overlapping of single-particle resonances may take place also in the Kondo experiments 
in QD-s [19], [501 EH [52], [53], where in order to increase the Kondo temperature the dot 
is usually sufficiently opened. Hopefully the unusual effects observed in some of these 
experiments may be also explained within our approach. 

Our mechanism of charging of the QD requires the existence of the broad level with 
T ^> A. The simple way to examine the relevance of our theory for the explanation of 
the experiment of ref. [1] will be to close the dot sufficiently in order to have T A for 
all levels. In this case the phase still increases by 7r at any resonance, but the correlation 
between peaks will disappear. (More precisely the pairs of peaks corresponding to adding 
of electrons with opposite spins onto one and the same level are still correlated, but the 
correlation between pairs should disappear.) Moreover, within our mechanism, a series 
of ~ (r/A) ln(f/cB/r) strong charging peaks in the conductance should have the same 
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height. This "coupling dependent" correlation of the peak heights seems also easy to 
measure. 

Next, we discussed the effects of spin. Here the most interesting predictions are the 
two-electron occupation switchings in the valley and sharp features inside the charging 
resonances. For the first time here we considered some fine structures of the population 
switching and their signatures in the signal of the charge detector. 

In the last section, the extra-strong sensitivity of the transmission phase through a 
quantum dot to the Kondo correlations was explained. The nontrivial behavior of the 
phase develops at temperatures large compared to T K and may be found analytically by 
means of a simple renormalization group calculation. The phase reflects the developing 
Kondo correlations in the regime where, although the spin of the impurity is not 
screened, the running Kondo coupling constant exceeds parametrically its bare value. 
New features of 0, which were not noticed in existing numerical simulations were found 
for the Anderson impurity model. 

Finally, we mention that the results of a recent experiment on the phase behavior in 
very small dots with not-too-large electron numbers [11] would fit very naturally with 
models that necessitate relatively large level-width to level-separation ratios for the 
universal phase correlations. This is because the level-separation obviously increases 
with decreasing dot sizes. In fact, for the small dots the universal behavior gives way to 
more random "mesoscopic" one. The universal to mesoscopic crossover was studied in 
Refs. [TT|[T8]. with emphasis on the relevance of Fano physics to this problem. Numerical 
studies of the interplay of interference, interactions and level-occupation switching have 
been presented in Ref . [60J . Systematic experimental studies invoking both the size and 
the opening of the dot as well as charge sensing together with transport, could shed 
more light on these fascinating phenomena. 
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